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ABSTRACT 

We develop a new method to estimate gravitational shear by adopting an 
elliptical weight function to measure background galaxy images. In doing so, 
we introduce a new concept of "zero plane" which is an imaginal source plane 
where shapes of all sources are perfect circles, and regard the intrinsic shear as 
the result of an imaginal lensing distortion. This makes the relation between the 
observed shear, the intrinsic shear and lensing distortion more simple and thus 
higher-order calculation more easy. The elliptical weight function allows us to 
measure the mutiplemoment of shape of background galaxies more precisely by 
weighting highly to brighter parts of image and moreover to reduce systematic 
error due to insufficient expansion of the weight function in the original approach 
of KSB. Point Spread Function(PSF) correction in E-HOLICs methods becomes 
more complicated than those in KSB methods. In this paper we studied isotropic 
PSF correction in detail. By adopting the lensing distortion as the ellipticity of 
the weight function, we are able to show that the shear estimation in E-HOLICs 
method reduces to solve a polynomial in the absolute magnitude of the distortion. 
We compare the systematic errors between our approach and KSB using STEP2 
simulation. It is confirmed that KSB method overestimate the input shear for 
images with large ellipticities, and E-HOLICs correctly estimate the input shear 
even for such images. Anisotropic PSF correction and analysis of real data will 
be presented in forthcoming paper. 
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1. Introduction 

Weak gravitation lensing analysis has been widely recognized as a very useful and unique 
method to study not only mass distribution of the universe, but also the cosmological pa- 
rameters(see for example Mellier 1999, Refregier 2003, van Waerbeke 2003, Munshi 2008). 
However the signal of weak lensing is small and thus accurate method of shear estimation 
must be required. Many method have been proposed so far and the most popular one is KSB 
method(Kaiser, Squires & Broadhurst 1995, Luppino & Kaiser 1997, Hoekstra et al 1998) 
which measures moments of the galactic light distribution and calculates change of moments 
by lensing distortion. The KSB method is successfully applied to many objects(Kaiser & 
Squires 1993, Broadhurst et al 2005, Umetsu & Broadhurst 2007, Okabe & Umetsu 2008 
and Okabe 2010) including cosmic shear which is distortions by large scale structure (Bacon 
et al 2000, Kaiser et al 2000, van Waerbeke et al 2001 and Hamana 2003). 

In particular the cosmic shear attracts much attention recently because of its potentiality 
to determine the so-called the cosmic equation of state, namely the relation between the 
energy density and pressure of the dark energy which will decide the destination of the 
universe. However the lensing signal of the cosmic shear is much smaller than the signal by 
cluster weak lensing and thus highly accurate measurement of shear is necessary in order to 
determine the cosmic equation of state in a percent level. There are several sources for noise 
in the estimation process of lensing distortion which we have to reduce as possible as we can. 
For example, sources have intrinsic ellipticities before affected by lensing distortion. Because 
we can not divide the estimated distortion into intrinsic ellipticity and lensing distortion, the 
intrinsic ellipticity behaves as a noise. We can reduce this noise by averaging over enough 
number of the estimated distortions because the intrinsic ellipticities are expected to have 
random orientations. In this way the noise from intrinsic ellipticity is controllable. However 
there are other noises which are uncontrollable such as noises coming from random photon 
count, smearing by atmosphere. We do not have methods to correct these effects completely, 
and insufficiency of correction for these effects becomes a systematic noise. Therefore we 
need to develop a method which can measure more accurately shape of background galaxies 
and correct more carefully the effect of the smearing. The precision and problems of KSB 
method have been studied (for example, Kuijken 1999, Erben et al. 2001, Hoekstra 2004, 
Hamana & Miyazaki 2008). Then new analysis methods (for example shap let (Refregier 2003, 
Refregier & Bacon 2003, Kuijken et al. 2006) and lensfit (Miller et al 2007 & Kiching et al 
2008)) and improvements of them (Recently, Viola et al. 2010, Melchior et al. 2010) have 
been developed. 

In order to test the accuracy of various shear measurement schemes, some testing pro- 
gram such STEP1 (Heymans et al 2005) and STEP2 (Massey et al 2006) have been developed. 



-3 - 



For example, the results of STEP2 showed that many of them reach the level of a percent 
accuracy. However, there are two major issues in this result. One is that nearly all methods 
need additional weights and extra corrections to achieve such an accuracy. Another is that 
most of the previous work are interested only in the averaged precisions, namely they are 
studied the averaged distortion of all images. Thus there is a possibility that high accuracy is 
the result of cancelation between the over-estimate and under-estimate of the distortion for 
different type of the images. Even if we have a set of lens images which has correct average 
value, variance increases(decreases) by overestimate (underestimate) of the shear. Especially 
a correlation between redshift and these over and/or under estimation affects seriously the 
cosmic shear estimation and the prediction for the cosmic equation of state. 

In this paper we develop a new method of the evaluating shear which is based on 
KSB method, but improves the above two issues. Our method is a natural development of 
our previous studies. One is oct-HOLICs (Okura et al 2007 and Okura et al 2008) which 
estimates higher order distortion "flexions" by using oct-pole components of HOLICs, and 
another is S2-HOLICs (Okura et al 2009) for precise measurement of shear by using high 
order spin-2 components and increasing information of image. We will generalize the HOLICs 
approach by adopting an elliptical weight, which we call elliptical weighted HOLICs(E- 
HOLICs) method. We will show that HOLICs method with an elliptical weight is able to 
measures the lensing distortion more accurately by weighting highly to brighter region of 
image than in the standard KSB method, and thus it can reduce effects of systematic error 
and random noise more effectively than KSB method. 

The organization of the paper is as follows. First we give a very brief review of weak 
lensing in section 2. Then in section 3, we introduce a new notations and definitions for 
HOLICs which make calculation a little easy, and then we introduce the concept of "zero 
plane" which is an imaginal source plane where all sources are perfect circle. The source 
plane is regarded as the result of an imaginal lens mapping from the zero plane. This allow 
us to separate the shear estimation into a process of evaluating ellipticity before smeared 
by PSF and a process of reducing noise of intrinsic ellipticities by averaging. KSB method 
is examined from our point of view in section 4. Then E-HOLICs method is presented in 
section 5. Tests of KSB method and E-HOLICs method using STEP2 simulation data are 
shown in section 6. Finally we summarize our result and give some discussion in section 7. 



2. Weak Lensing convergence, shear and distortion 

A gravitational deflection of light rays can be described by the lens equation as 

P = 6-d*{j{6), 



(1) 
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where 6 = 61+162 and (3 = Pi+ifa are angular positions of the image and source, respectively, 
d = di + i&2 is a complex gradient operator that transforms as a vector, so d' = de 1 ^ with 
<f> being the angle of rotation, and ip(6) is the effective lensing potential, which is defined 
by the two-dimensional Poisson equation as V 2 ip(6) = 2n(6), with the lensing convergence. 
Here the convergence k = S m S~J t is the dimensionless surface mass density projected on 
the sky, normalized with respect to the critical surface mass density of gravitational lensing 
S C rit = (c 2 D s )/ (4nGD d D ds ), where D d , D s , and D ds are the angular diameter distances from 
the observer to the deflector, from the observer to the source, and from the deflector to the 
source, respectively. 

The lensing convergence k is expressed as 

k = ±dd*i>, (2) 

where * denotes a complex conjugate. Similarly, the complex gravitational shear is defined 
as 

7 = 7! + ry 2 = -ddi) (3) 

which has a property of spin-2. Note that a quantity is said to have spin-s if it has the same 
value after rotation by 2ir/s. 

Because these convergence and shear have a relation in Fourier space as 

«(*) = ®~$~% klk2 m, (4) 

therefore mass distribution of lens object can be obtained from shear distribution. 

Here we define a reduced shear g and a complex distortion 8 which appear in following 
sections as 



3. New Notations, New Definitions and Zero plane 

In this section, first we introduce complex forms in displacements from centroid of image 
and moments of image. These complex form makes equations much simpler and higher-order 
calculations easier. Next, we introduce "Zero image" which is imaginary circular image and 
becomes "source" by an imaginary distortion. In estimating shear, intrinsic ellipticities of 
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background galaxies complicate equations, because an ellipticity of distorted image is not 
a simple sum of the intrinsic ellipticity and distortion. By treating the intrinsic ellipticity 
as an imaginary distortion from the zero image, the relation between the observed shear, 
the intrinsic shear and the lensing distortion becomes a simple form(see eq.(31) below), and 
we are able to divide the shear estimation process into the estimation of individual distor- 
tions(which contains intrinsic ellipticity) and the reduction of noise of intrinsic ellipticities 
by averaging. 



3.1. Complex Notation 

We use complex notation for displacements and moments of image which are expressed 
by subscripts in KSB method. 



3.1.1. Complex Displacement 
Let 9 = 9i+ i9 2 be a centroid of image and a displacement from this point is notated as 

A9 = 9-9 = (9 1 -9 1 ) + i(9 2 - 9 2 ). (7) 
Then a combination of displacements is written as 

»r N+M N—M 

A9^ = A9^-A9* — , (8) 
where " N" means order of length and " M" means a spin number. 



3.1.2. Complex Moments and HOLICs 
Complex moments of an image without weight function is defined as follows. 

ZZ{I) = J d*9I(9)A9 N M , (9) 

where 1(9) is the brightness distribution. Similarly the complex moments with weight func- 
tion W which has an ellipticity value same as d is defined as 

z-d, g) S j ttwuizw , (10) 



-6 - 



where "N" and "M" mean order of length and spin number respectively. Because distortion 
makes ellipticity, distortion and ellipticity have same meaning in E-HOLICs method each 
other. For example, comparing moments with weight function with KSB notations becomes 

Zl(I,0) = Qn-Q 2 2 + 2iQ 12 (11) 
Z 2 (I,0) = Q n + Q 22 , (12) 

where Qij is the components of quadrupole moment in KSB notation. 

In addition, we introduce HOLICs s"H^ which is the complex moment Zj^(I) normalized 
by Zx(I) as follow and notate as 

<(/,^§|g (13) 
Similarly we introduce elliptically weighted HOLICs(E-HOLICs) as 

nUi,z L K ,s)^^§-. (14) 

For notational simplicity, we write a combination of HOLICs as 

nm, Z L K) 8) 



njN I njP 
rt M -+- ft 



O 



n/R 



(I, Z L K , 8) 



(15) 



3.2. Source, Image and Zero image without weights 

Here, we define an imaginal brightness distribution "Zero image" which allows us to 
calculate distortion easily, and present the relations between "Zero image", "Source" and 
" Image" . 



3.2.1. Source and Image 

First, we present a relation between Source and Image which can be written by conver- 
gence and shear. 

A relation between A/3 and A9 which are complex displacements in source plane and 
image plane, respectively, can be written in matrix expression as follows, 

( Z\ ) - ( ^ A ) ( £ ) - <-> ( ';f A ) ( £ ) • 
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where g L = 7 L /(1 — k) is reduced shear, and this relation can be written in our complex 
notation easily as 

A/3 = (1 - k)A9 - 7 L Ar = (1 - «) (A6> - g L A9*^j . (17) 

By writing brightness distributions of Source and Image as I s '(/?) and I L (9), respectively, 
we obtain the following relation for complex moments of source and image. 

35(/ fl ) ss J rf 2 /3/ 5 (/3)A/3^ = (1 - «)" +2 (l - |/| 2 ) | <Pei L (e)(M - g L A9*) N M . (is) 



Source and Zero image 

Next, we present a relation between Source and imaginal Zero image , where I z and A/3 
are a brightness distribution and a complex displacement of a zero image, respectively. The 
zero image is defined as a brightness distribution which has no non-spin-0 moments, and we 
can relate zero plane and source plane by imaginal shear (= g 1 , there are no imaginal 
convergence) as 

A/3 = A/3 - 7 J A/5*. (19) 

Here we assumed that there is no imaginal convergence, so the imaginal reduced shear is 
defined as 

g 1 = i 1 . (20) 

From the above definitions, the complex moment Z|(/ z ) of a zero image can be calcu- 
lated as 



= Z 2 2 (I z ) = J d^I z ((3)Afe (21) 

= (1 - IsT) / d^I s {P) {API - 2g J API + g»Ap?) (22) 
= (1 - 1ST) (^(/ 5 )-2^ 2 (/ 5 ) + |/| 2 Z 2 2 (/ 5 )) (23) 



= (1 - 1ST) (1 + i*T) (^ 2 2 (/ 5 ) - r^7p^o(^)) 



(24) 



= (1 - \ 9 r) (Z 2 2 (I S ) - S'ZZil 3 )) ■ (25) 

where we have used the fact that g 12 Z\*{I s ) = |g 7 | 2 Zf (J 5 ), namely Z\{I S ) and g 1 have 
same phase angle because Z\{I S } is distorted by g 1 from a zero image. Then we can obtain 
ellipticity easily as, 

Hi(I s ,Z 2 ) = 5 I . (26) 
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Therefore we can simply use 5 1 and g 1 instead of intrinsic ellipticity %|(/ 5 , Zq). 



3.2.3. Image and Zero image, and Lensing shear 

From the above relations, Image can be obtained by two distortions by g 1 and g L from 
Zero image. Complex displacement in image plane is distorted as 

A/3 = (A/3 - /A/3*) = (1 - k) (AB - g L A9* - # Z (A0* - g L *A9)) (27) 

= (1 - K ) (1 - /V) (A0 - g c A9*) , (28) 

where 

f - (29) 

1-gIgL* V J 

is the combined shear. From similar calculation of eq. fT2~T|) . we obtain 



Ht(I L ,Zl) = 6 c = T ^ Jr2 . (30) 



Because g 1 has random phase angle, Lensing shear g L can be obtained as shear which satisfies 
following equation 

'=<^(iW>' (31) 

This is our basic relation. 

By introducing the Zero plane, we can divide lensing analysis into two steps which are 
the evaluation of g° and then the reduction of noise in intrinsic shear by eq. (!3T|) . Because 
we don't use any assumption and the property of E-HOLICs for introducing eq. (l3~T|) . this 
equation can be used for any methods of weak lensing shear estimating. Therefore our 
purpose is to obtain g c in this paper. And here and after, g = g c and S = 8 C for simply 
notation. 



3.2.4- Shear Distortion for Higher order moments 

We demonstrated that the effects of lensing distortion for order-2 spin-2 moment is 
7-i\{I L , Zq) = 6. From similar calculation of eq. f)30p . relations between distortion and higher 
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order HOLICs are obtained as 





= s 


(32) 


Zq) 


3 . 
= 2 + P* 


(33) 


^SJ-, Zq) 


3 A 2 


(34) 




4 + <5 2 


(35) 


2 + 35 2 


Ut(I, Zq) 


5 S 2 

2 + 3d 2 


(36) 


nt{i,zt) 


= 5 5* 
2 + 35 2 ' 


(37) 



where 5 = \8\. These equations will be used in forthcoming papers. 



4. KSB method using zero plane in new notation 

In this section we review the standard KSB method using zero plane and in our new 
notation. For simplicity we suppose 6 = (e.g. 6^ = A0 M ). The detailed explanation and 
calculation can be seen in Bartelmann and Schneider 2001. 



4.1. Weight function 

In our approach it is natural to assume a circular weight function in the zero plane rather 
than in the image plane as in the original KSB method. This will give us the unaveraged 
expression for PSF corrected shear (see equation. (184]) ) and the lensing shear will be obtained 
by averaging using equation (I3T1) . The relation between weighted complex moments and 
complex distortion in our case becomes as 



Z»(I*,0) = I #pi z {fi)p»W ( ^ ] (38) 



k ) N+2 (1 - \g\ 2 ) f d 2 9I L (9) (9 - g9% W ( ° l R *} 5 *° 1] ) . (39) 
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By calculating complex moment Zf, a relation between ellipticity and distortion can be 
obtained as 

= Z 2 (I z ,0) = [d 2 pi z 0)%W ( |f ] (40) 



\g\ 2 ) J d 2 6I L {0) (6 - gP)l W ^gM ) (41) 

= (1 - - \g\ 2 ) (Zl(I L , S) - 2gZ 2 (I L , 6) + g 2 Z 2 2 *(I L , 6)) (42) 
= (1 - - \g\ 2 ){\ + \g\ 2 ) (Z 2 (I L , 6) - SZ 2 (I L , 6)) (43) 

« (1 - k) 4 (1 - M 4 ) [z\{l\ 0) - £- 2 Z't(I L , 0) - ^1(I £ , 0) - SZ 2 (I L , 0)) ,(44) 

where we have expanded the elliptical weight function and kept only 1st order in 5 because 
KSB method use circular weight function. 



M '"W(%), (45) 



2<7 2 \<ji 

and Z' is the complex moment with weight function dW(x) / dx instead oiW{x). Therefore, 

zl o) « (1 + ^n'i(i L , zl o)) s + ^ o)) <r (46) 

_ C KSB^ S + C KSB^ 5 * ) (47) 



where is defined as 



0)" 

A complex distortion is thus calculated as follows. 

sKsb ^|(/ L , Zl 0)C (I L , 0) - Hf (I\ Z 2 , 0)C A (I L } 0) 
\C {IL,Q)\ 2 -\C*{I L ,0)\ 3 



(49) 



If we consider only 1st order of combinations in C and H^, we can neglect this terms like 
H 2 *C±. Therefore we can obtain shear simply as 

eKSB _ 7-Ll{I L ,Zlti) 

6 ~ cJFJ) (50) 



Because the approximation in eq. ()49p is effective under the assumption of small S, 
the equation has systematic error in high 6. For example, if image and weight function 
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are Gaussian with same size and have ellipticity 8 tr ue and 0, respectively, we can calculate 
complex moments as 

Z» = J d 2 66 N M e^l° 2 e- e «l° 2 (51) 

and obtain 

nl = (52) 

n KSB 8 — hb\ Tue , 
°° ~ 4(4 (53) 



2 



n KSB _ Sdtrv r , - . , 

° 4 ~ 4(4-^)' (M) 
Therefore the estimated distortion S estimate is derived as 

1 - ft 2 I A 

s s x u true I ^ 

"estimate utrue -. ^2 • V / 

^true 

Because KSB method corrects only the numerator of H%, one can over estimate the distortion 
in general. 



4.2. PSF correction 

Because the observed images are smeared by atmospheric turbulence and so on, we must 
correct these effect to obtain true distortions. This effect can be expressed by Point Spread 
Function(PSF) and can be divided into isotropic PSF which reduce ellipticity of lensed image 
and anisotropic PSF which add a extra ellipticity. KSB method corrects these two effects 
using moments of star images which is smeared by same PSF. 

An observed image I obs (9) is smeared by PSF P(9) from lensed image I L (f3) and ex- 
pressed as 

I obs (9) = J dpi L {P)P{B-P). (56) 

KSB assumes that PSF is expressed as a convolution of an isotropic part P lso and an 
anisotropic part of q as follows. 



P{Q) = / d 2 i)P lso {4))q{B - V). 



(57) 
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The moments of anisotropic PSF is assumed to satisfy the following relations. 

4 = 1, 

q\ = o, 



qg = 0. 



where moments of anisotropic PSF are defined as 



d 2 ee N M q{6) = q N M , 



Let P so be an image which is smeared by only isotropic PSF as 

P so {B) = J d 2 /3I L ((3)P iso (6 - (3), 

then eq. fl56l) becomes 

P hs (0) = [ d 2 i)I iso {i))q{9 -if)). 



(58) 
(59) 
(60) 

(61) 

(62) 
(63) 



4-2.1. Anisotropic PSF correction 

Here, we explain a method of anisotropic PSF correction. 
Moments of I obs are described by integral of P so and q as 

d 2 9f(9)I obs (6) = [ d 2 6 J d 2 i)f(d)I iso (i))q(d -if>) 

d 2 <j> [ d 2 ^f{<j> + ^)I is0 ^)q{<j>) 
'9 2 -Re[d*e 2 } 



e N M w 



Expanding f((f) + ip) up to 1st order of 5, Z 2 (I , 0) is obtained as 



Z 2 (I obs ,0) 



Z 2 2 +ql[Z [ 



a 



2 „,2 1 7 „4 



2* 7 //4 
Q2 L 4 



(P so ,0), 



Then we have an approximated expression for 7i 2 (I ,0). 



U\{l ob \ Z 2 ,0) 



%2 + 92 ( ^0 + T- T^J 



1 «i ) + <zr< 



^(/ iso , Z 2 , 0) + q 2 2 P KSB (I obs ) + ql*P* SB (I obs ). 



(64) 
(65) 
(66) 

(67) 



(P s °,Z£,0) (68) 



(P bs ,Z 2 ,0)(69) 
(70) 
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Similar calculation is applied for star image I* and we obtain 

n 2 2 (r obs , z 2 , o) = q 2 2 p* CSB (i* obs ) + ql*pf SB (r obs ) (71) 

This will give us the following expression for q\. 

2 _ u\{r oh \ z 2 , o)p£ CSB (i* obs ) - nf{i* obs ) z 2 , o)pf SB (i* obs ) 

If we use 1st order of the above equation, V? 2 P^ SB may be neglected, and we have 

2 _ Hl(P obs ,ZlO) 

^2 pKSB(j*obs^ ■ v'-Jj 

where P^ SB is trP which is often used instead of matrix P in weak lensing analysis with 
KSB method. 



4-2.2. Isotropic PSF correction 

Next, we explain the method of isotropic PSF correction which reduce ellipticity of 
images. Let I(/3) be a brightness distribution in zero plane. According to the conservation 
of brightness, this coincides with P so {6) in image plane, and we can describe as 



r»>(0) = m = J ww)P{p - v). (74) 

Because P is not isotropic function, we write P in the following way. 

P0) = J d 2 4>P is °$)q0-4>). (75) 

Let P so be a brightness distribution which is smeared by only P tso (ip), and 10) can be 
described as 

10) = J d 2 ^P so 0)q0-^). (76) 

Because the anisotropic part q is made by shear, we can assume that q^ and q\ are dominant 
in q. Therefore we can correct anisotropic part similarly as above, so we obtain 

= H 2 2 (P s °,Z 2 ,0) » n 2 2 (I,Z 2 ,0)-P KSB (I)q 2 2 -P, KSB (I)qf (77) 
« n 2 2 (I, Z 2 , 0) - P^ SB {I obs )q 2 2 - P« SB (I obs )q 2 2 *. (78) 
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We need the relation between 7-LI(I,Zq,0) and H%(P so , Z$ , 0) which is obtained by similar 
calculation as eq.f )46|) . 



H 2 2 (I,ZlO) « Hl{f so ,Zl^)-C^ SB {f so )8-Cf SB {r o )8* (79) 
« 7*1(1*°, Z 2 , 0) - C^ SB (I obs )8 - Cf SB (I ohs )8*. (80) 

Thus 7/|(/* so , Zq, 0) can be described by 8 and q\ as follows. 

^(I iS °, Z 2 , 0) « fifB (/ o65 )« + Cf SB(jobs) S * + pKSB ( p6 S) ~2 + plCSB(p^~2* ^ 

g 2 , can be obtained from star images as 

= ^ (/ «so ; Z 2 ; Q) w C * SB (r)8 + Cf SB {T)8* + Pf SB (P)q 2 2 + Pf SB {P)qf (82) 

[pKSB pKSB _ QKSB*pKSB^ § _|_ ^jKSBpKSB _ QKSBpKSB^ 

(J*), (83) 
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Substituting this expression to eq. f[8"Tj) . we can estimate complex distortion. Especially in 
1st order, we can use the following simple equation to estimate shear. 

£ ~ fiKSB = Hl(P so , Zl, 0) 

~ rKSB(jobs\ r KSB(T*\ PZ [SB {I° bs ) ' K ' 

°o K 1 ) ~ °o K 1 ) pits By*) 



5. Elliptical weighted HOLICs 

In this section, we demonstrate an improvement for the overestimation of shear in KSB 
by adopting an elliptical weight function. In this paper we will not consider anisotropic PSF 
correction which will be treated in the forthcoming paper. 



5.1. Motivation and Merits of using E-HOLICs 

The one of the major concern in weak lensing is to measure the cosmic equation of 
state by cosmic shear. It is needless to say that this will require us to have highly accurate 
method to evaluate the gravitational shear. Although KSB method is very powerful and 
convenient, it uses only the first-order expansion of the weight function as we demonstrated 
in the previous sections, and thus does not have enough accuracy in such a purpose. For 
example, KSB has 7.5% bias if the lensed object has ellipticity 0.3 which is typical value 
of intrinsic ellipticities. Moreover, the use of a circular weight for elliptical images in KSB 
method is not an efficient and accurate measurement for galaxy shape. 
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Shear Testing Programme 2(STEP2) is a project for testing precision of weak lensing 
measurement (See Massey et al. 2007). It uses several distortions and 6 PSF patterns for 
same source objects, and tested precision of various shear measurement methods including 
KSB. Figure 6 of the paper by Massey et. al (2007) shows that the precisions are about 
1 percent in most of the methods they tested. However, Table 4 of the same paper shows 
that artificial corrections are needed to obtain such a precision. It is natural to imagine 
that such artificial corrections will change by every real data and every shot of observation, 
because the artificial corrections corrects the systematic error which are due to size, MAG 
and SN of objects, size of PSF and so on. We need some sort of correction in each data, but 
it will be hard to have universal corrections independent of various observational conditions 
and various PSFs. For reducing such hardness, it is thus desirable to develop a new shear 
measurement methods which need fewer or/and less artificial corrections. 

One of the reason to need such corrections is the use of insufficient expansion of the 
weight function as demonstrated in section 4. This lead us to adopt the elliptical weight 
function to measure the galaxy shape with the distortion given by the lensing distortion. 
This will give us a relatively simple relation between the observed, intrinsic shear and lensing 
distortion. In addition, because using same ellipticity of image for weight function increases 
SN of measurement, we can reduce random count noise. 

In this section we show the detail of our approach to estimate a lensing shear by using 
a elliptical weight function. 



E-HOLICs method is defined that the weight functions to measure the shape must have 
no ellipticity in zero planes. This directly gives us the following simple expression for spin- 
2 E-HOLICs without any approximations according to the similar calculations leading to 
eq.flU). 



5.2. Distortion with weight function 



Zj(I L ,S) 
Z 2 (I L ,S) 



H 2 2 (I L ,Z 2 ,d) = 5. 
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Moreover, the relations between complex distortion and higher order E-HOLICs can be 
obtained similarly as eq. fl32|) . 

U 2 2 (I,Z 2 ,8) = 6 (86) 

= (87) 

WftfX*) = ^-2 5 " (88) 

W^X*) = (89) 

«2('X*) = ( 9 °) 

wS(/X*) = ^f2 63 ( 91 ) 



5.3. One dimensional form 

In next section, we demonstrate how to correct isotropic PSF. However because our 
weight function has ellipticity, PSF corrections are more complicated than that of KSB 
methods. Therefore it is essential to reduce the equations to one dimensional form in order 
to make the method tractable in real observations. We will explain what is meant by one- 
dimension in the below. 

First, we divide complex distortion and complex moments into absolute value and phase 
angle. Let <fi$ be a phase angle of 5, so 

8 = 5e i(t>s . (92) 

And <t>% be a phase angle of H^, so 

Ul{l\ Zl 6) = H»{I L , Zl (93) 
where = \7i^\. Owing to eq. fl85|) . we can obtain 

H*(I L X,8)=5, (94) 
02 = 05- (95) 

Next, we calculate the correlations between phase angle of complex distortion and com- 
plex moments using STEP2(we describe STEP2 simulation images in following section). The 
figures [Q, [2] show the difference of the phase angle between higher-order(N=4, 6) sipn-2 E- 
HOLICs and complex distortion as a function of ellipticity of images. These figures clearly 



-17- 



show a high correlation (02 — (ps) between them independent of ellipticity of images. On the 
other hand, figure [3] shows the difference of phase between spin-4 HOLICs with N=4 and the 
complex distortion as a function of ellipticity of images. This shows also a high correlation 
(204 ~ <f>g) at least for images with high ellipticities. Fortunately, T-L\ appears always as a 
spin-2 combination with other complex distortion (for example S*TLf), so terms with b*li\ 
are always 2nd order or more higher order. Therefore, although the phase of H\ has a low 
correlation with that of the lensing distortion if object have low ellipticities, it can safely 
be neglected in our analysis as far as we are interested only in the lowest order calculation. 
Similarly terms containing other higher-order moments with higher spin are always smaller 
than 2nd order. Thus we employ the following approximation for the phase angles of various 
E-HOLICs. 

Ki ~ y</>2 = Y 0«5- (96) 

This approximation allows us to reduce complex equations for the lensing distortion into real 
equations of its absolute value. 



5.4. Isotropic PSF correction 

E-HOLICs method uses similar idea of KSB for PSF correction. In this paper we 
present only isotropic PSF correction, and we will present the anisotropic PSF correction by 
forthcoming paper. 

E-HOLICs method requires Gaussian weight function in the PSF correction, so 

°0 ~ Re WV \ -(flg-Re[a*0i}])/ff 



W 21 = e -^5-^L° W. (97) 



First, we suppose that an isotropic PSF P lso (f3) in the zero plane can be divided into isotropic 
part and anisotropic part in the image plane as follows. 

P iso 0) = P{6) = J d 2 <j)P iso (<P)q(9 - 0). (98) 

where P(0) has the ellipticity which has same value of complex distortion, and P lso {9) is the 
isotropic part of real PSF. 

Let P so be the brightness distribution which is smeared by P lso (f3), then the complex 
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moments of P so can be expressed by the complex moments of P so and q as 

J d 2 p f0)P*°0) = Jd 2 Pf0)Jd 2 ^P iso 0-^)I z ^) (99) 
= (1 - k)\1 - \g\ 2 ) J d 2 9f(9 - g9*) J d<pI L (<p)P(9 - <p) (100) 
= (1 - nf{l - \g\ 2 ) J d 2 9f(9 - g9*) J d V I L { V ) J d 2 <t>P is °{9 -if- 0)g(0)(lOl) 
= (1 - k) 2 (1 - \g\ 2 ) J d 2 i> J d 2 <j>f((l - - #* + <f> - g4>*))P so {m{4>im 
where /((l — k) (i/j — gip* + <p — g<p*)) is 

/((l - k) - #* + <f> - g<f>*)) (103) 

g (1 _ Kf{ , _ ^ + , _ g , ruW ( (* + - r : + 0)1] \ ; (104) 



and W is expanded by using the property of Gaussian as 



tr[ ty + fl§-Re[J*ty + flj] > | (1()r)) 



i - 1 (V + r<P - w - <w - ^vi - ^r) (los) 
f^ * + " 5 * n - sr<t) * " ~ \ 54>2 ) ^ ■ (109) 



Here, the complex moments of q is defined as 



f d?M»W ($) q(0) 



where # M satisfies g{J = 1 and qff = 0(7V 7^ 0), and we assume that and are dominant 
and we ignore other moments of q. Moreover, because the phase angle of e 1 ^" (defined as 
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q\ = q 2 e l< ^ q ) is same as 0<j, eg. (199]) reduces to one dimensional form and we obtain 



Zi(P so ,0)^(l-Kni-\g\ 



Z 2 ~ $ Z 



111) 



+ Ql (zl - ± ((4 + 36 2 )Z 2 - 76* Zl) + ^ ((1 + 25 2 )Z 4 - 6*(2 + 5 2 )Z$ - 8Zf + 6* 2 2$})j) 

(I lso ,d) (113) 

Hi- 5 + q 2 (tf ° - 1 ((2 + 35 2 )H 2 - 55tf 2 ) (114) 



$ ( -^( 3 ^o - ^Z 2 2 ) + ^(3<5% 4 - 3<5Z 2 4 + d 3 Zt + Z 4 4 



:i-K)\i-\g\*)z 2 (r s °,dy< 



((1 + h62)H » ~ 26(3 + 62)H2 + + 62)Ht) ) 



(P s ° } zl 8) 

= (1 - k)\\ - \gf)Z 2 (P so , *) e «*» [H 2 2 (r°, Zl 6) -6 + q 2 2 P E (P so , S) ^ 

where P E is a combination of absolute values of HOLICs and d, thus P E is real. Therefore, 
from the above equation, we have the equation for 5. 

5 » H 2 (P S0 , Zl S) + q 2 2 P E (P so , 6). (117) 

This equation can be apply to star objects to have an expression for q 2 . 

S -H 2 (P is °,Zl5) 



(115) 
(116) 



?1 
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P E (I* iao , 6) 

where P lso is circle but H 2 does not vanish because our weight function has ellipticity. 
Finally, the lensing distortion can be evaluated by solving the following equation 

P E (I iso , S) 



5^5 E ^H 2 (P so ,ZlS) 
The complex distortion then is obtained as 



l2 I TISO r/2 



pE^J*isc 


,*) 


as 




P E {I iso , 


*) 




,6) 


P E (P S0 , 


S) 



P E (P is °,8) 



{5-H 2 2 (P^,Z 2 ,d)) 

(5-H 2 2 (P* s °,Zl5)))e> 
(S-H 2 2 (P iso ,Zld)) 



(119) 

(120) 
(121) 



5.5. Ellipticity of weight function 



Eq. (11191) requires us to use the ellipticity of lensing distortion. However an observed 
object has a different ellipticity from the lensed ellipticity due to PSF. Therefore we can not 
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obtain lensing distortion directly by observation, and thus we use temporally the observed 
ellipticity as a natural first step. Then there are two ways to estimate the lensing distortion. 
One is iteration. Namely we use the observed ellipticity for elliptical weight and re-estimating 
lensed ellipticity, and iterate measurement by using re-estimated ellipticity until the result 
converges. Otherwise we could expand the weight function for all moments of Eq. (I119j) as 

W (M) (122) 
' {(S - 8w? 0t + (S- 6 W ) (<T - S* w ) # 4 + (6* - 8* w ) 2 6t) ) , (124) 



4a 4 



where 8w is an ellipticity we use for weight function. It seems reasonable to think that 
coefficients of this expansion are smaller than these of KSB method. By using observed 
ellipticity for $w, systematic error due to this expansion will be smaller than that in KSB. 

We have examined both method and found the iteration needs more time than the 
expansion. Thus the above expansion of the weight function may be more useful in real 
observation because we must estimate PSF effects for many star objects. 

Two dimensional fitting set of with elliptical weight functions which have all phase 
angle and ellipticity or using only some star objects which surround each distorted object 
which we try to correct PSF effects will be effective in saving time to calculate. 



6. STEP2 simulation test 

In this section, we perform tests of our methods and compare our results with KSB 
method by using STEP2 simulation data. Because we presented only isotropic PSF correc- 
tion, we use data set of psf-F which have no anisotropic part. There are 64 data set of STEP2 
simulation data which are distorted by different distortion respectively and have a rotated 
image. By using these rotated data, an effect of intrinsic ellipticity can be canceled (but a 
systematic error of intrinsic ellipticity can not be canceled completely). In our evaluation 
of shear we have ignored denominator of eq. (l3T|) when averaging for neglecting systematic 
error. In real analysis, we should use iteration to consider higher order in the equation. A 
perl scripts (Umetsu Keiichi private communication) and IMCAT were used for detecting 
objects. In averaging with the normal and the rotated data, the number of matched objects 
in No. 36 field in STEP2 is about 1/10 of the number of matched objects in other fields(may 
be due to some strangeness in detection), thus we neglected the result of the field. 
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6.1. Star selection 

We selected star objects by the condition 1.9< <2.2(pix), 20<MAG<23. Generally, 
two dimensional fitting is used for PSF parameters in real weak lensing analysis, because PSF 
effect slightly changes in a data field. However, constant PSF is used in STEP2 simulation 
data, so we can correct PSF effects by using only one star object. Here we use one of good 
star which has SN^500 for PSF correction. In order to see the detailed behavior of the 
estimation depending on the choice of the objects, we classify the objects into the following 
classes. 



6.2. good objects 

First we define a class which contains only good objects. Good objects is defined to 
be objects which are large (r^ >3.3(pix)) and have high SN (SN>20) and small intrinsic 
ellipticities( | intrinsic | <0.35). We do not use any limit on MAG(e.g. 20<MAG<26). These 
conditions reduces objects used for the shear measurement where KSB method have a number 
density 1.50(/min 2 ) and E-HOLICs method have a number density 1.68(/min 2 ). Large 
images have fewer effect from PSF and high SN images have fewer effect from random 
count noise. The effect of expansion of weight function will not cause large error in shape 
measurement for images with small intrinsic images. Thus the estimated shear using only 
good object is expected to be close to the input value. 

Figure H] shows results of shear estimation by KSB method and E-HOLICs method using 
only good objects. In these figures the horizontal axis means inputted(answer) shear and 
the vertical axis means a difference between estimated shear and inputted shear. 

The result is described by two parameters m and c defined as 7 e = m^i + c where 7 e is 
the estimated shear and 7, is the input shear. KSB gives m=0. 03868 and c=0. 000232, and 
E-HOLICs gives m=— 0.0242 and c=0. 000083. Both give satisfactory result even if we used 
only a few number of images. Comparing with both methods, we can see KSB method gives 
slight over estimation for objects with large ellipticities. 

6.3. intrinsic ellipticity, SN and size dependence 

Next we would like to see the effect of intrinsic ellipticity, SN and size of objects sepa- 
rately in the shear estimation. 

First we relax the condition for intrinsic ellipticity keeping other conditions. We include 
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object with \5i n t r insic\ <1- This increases the number of objects up to 2.48(/min 2 ) for KSB 
and 3.41(/min 2 ) for E-HOLICs. Figure [5] shows the results where KSB gives m=0.3108 and 
c=0.000737 and E-HOLICs gives m=0.0044 and c=-0.000104. We can see that KSB method 
overestimates the input value strongly, The estimation by E-HOLICs method does improve, 
but the degree of improvement in this test is not explained by the increase of number density. 
It seems the this increase is explained by statistical variance. This is an evidence of that the 
expansion of weight function in KSB method is insufficient. Moreover, the result by KSB 
method has larger error bar compared to E-HOLICs method. This is because systematic 
error due to ellipticity does not cancel completely. 

To see the effect of size dependence in the shear estimation, we relax the condition 
of the size keeping other conditions fixed. We include objects with 2.2< on the top of 
good objects, and this increases number density of the objects up to 3.85(/min 2 ) for KSB 
and 4.05(/min 2 ) for E-HOLICs. Figure |6] shows results where KSB gives m=0.0277 and 
c=0.000231 and E-HOLICs gives m = -0.0075 and c = -0.000233. We can see there is no 
strong effects from dependence in both method. The reason for huge improvement may 
also be due to statistical variance. 

Finally we examine the effect of SN in the shear measurement. We therefore relax the 
condition for SN keeping other conditions fixed. We include objects with SN>2 on the top of 
good objects, which increases the number density up to 4.99(/min 2 ) for KSB and 5.38(/min 2 ) 
for E-HOLICs. 

Figure [7] shows results where KSB gives m=— 0.0539 and c=0. 000105 and E-HOLICs 
gives m=— 0.0604 and c=— 0.000553. We can see both methods underestimate the input 
value slightly and gain larger errors. This may be explained by the fact that random count 
noise reduce the observed ellipticities by adding random ellipticities. Because the degree of 
reduction by KSB method is larger than that of E-HOLICs method, one can conclude that 
KSB method is more affected from random count noise than E-HOLICs methods. 

We also found that the use of stars with low SN underestimates the input value, because 
random count noise disturb a precise measurement of an effect of PSF. Unfortunately, using 
only high SN stars is incompatible with small scale PSF correction. Because this problem 
is related with anisotropic PSF correction, we will address this problem in the forthcoming 
paper. 
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6.4. ALL object 

Finally, figure [8] shows the test using all object without any conditions mentioned above. 
Namely we use objects with 2.2< r^A|5| <1, 4<SN. This will increase the number density 
available for the measurement up to 23.70(/min 2 ) for KSB and 24.42(/min 2 ) for E-HOLICs. 
The result shows that m=0.0099 and c=-0.000261 in KSB method and m=-0.0014 and 
c=— 0.000863 in E-HOLICs method. We can see both methods give good estimation for input 
shear in this case. However, it should be mentioned that the accurate estimation obtained 
by KSB method is the result of fortunate cancelation between overestimate for images with 
large ellipticities and underestimate for images with low SN. On the other hand, because 
E-HOLICs method improves only overestimation due to the expansion of weight function, 
E-HOLICs method has a tendency to slightly underestimate the input shear. 

7. Summary 

We have developed a new method to estimate gravitational shear, called "E-HOLICs 
method" based on KSB method by introducing the elliptical weight to define multipole 
moments of galaxy light distribution. In E-HOLICs method we use the lensed ellipticity 
for the ellipticity of weight function, and thus it avoids a systematic error coming from an 
expansion of weight function which is done in previous approaches for shape measurement 
including KSB method. Furthermore the elliptical weight measures the galactic shape more 
effectively than with circular weight function, and it is expected to reduces random count 
noise. 

We have also developed the isotropic PSF correction for E-HOLICs method and were 
able to show that the equation which governs the correction reduces to real polynomial for 
the magnitude of the lensing distortion. This is possible because there is a high correlations 
between the phase angle of various order of E-HOLICs with same spin. For the application to 
real observation we need to correct anisotropic part of the effect PSF which will be discussed 
in the forthcoming paper. 

We have also introduced a new concept of "zero plane" in weak lensing analysis which is 
an imaginal source plane where shape of all sources is perfect circle, and the intrinsic shear is 
regarded as the result of an imaginal lensing distortion. This makes the relation between the 
observed shear, the intrinsic shear and lensing distortion simple, and higher-order analysis 
easy. The idea of zero plane may be used in any shear estimation scheme. 

We have performed various tests to show the improvements of shear measurement by E- 
HOLICs method using STEP2 simulation data, and compared the results with those obtained 
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by using KSB method. We have examine effects of ellipticity, size and SN of images separately 
in both method, and found that E-HOLICs method can reduce the errors coming from these 
effects. In particular, we have confirmed that KSB method overestimates the shear if images 
with large ellipticities are used in the analysis, and E-HOLICs method does not have such 
tendency. 

Although E-HOLICs method has a potentiality to accurately measure the shear, it 
requires more time for the estimation of the shear than KSB method in general, because 
it requires iteration to measure the multipole moments of background galaxies. Therefore 
it will be necessary to develop more effective method for the shape measurement when we 
apply E-HOLICs method to planned large scale cosmic shear observations. 

We thank K. Umetsu for a useful discussion and providing his scripts, T. Okamura for 
helping us in calculation, T. Hamana and M. Oguri for a useful discussion. This work is 
supported in part by the GCOE Program "Weaving Science Web beyond Particle-matter Hi- 
erarchy" at Tohoku University and by a Grant-in-Aid for Scientific Research from JSPS(Nos. 
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Fig. 1. — Plots of a ellipticity versus a difference of phase angles between V\ and %f of 
STEP2 objects. The horizontal axis means an absolute value of ellipticity, and the vertical 
axis means a phase angle (rad). Objects that of the difference is close to have high 
correlation in phase angles. 




Fig. 2. — Same plots as Figure [T]except the difference. The difference of phase angle between 
and 1-L\ is shown in this figure. 




Fig. 3. — Same plots as Figure [T]except the difference. The difference of phase angle between 
%\ and twice V\ is shown in this figure. 




Fig. 4. — Plots of results of test of 63 fields of STEP2, The left side shows a result of 
KSB method and the right side shows a result of E-HOLICs method. The horizontal axis 
means value of inputted shear (answer), and the vertical axis means a value of difference 
between estimated shear and inputted shear. The subscripts 1 and 2 mean components of 
directions. KSB method has a number density 1.50(/min 2 ) and obtained fitting parame- 
ter 0.03867+0.000232. E-HOLICs method has a number density 1.68(/min 2 ) and obtained 
fitting parameter -0.02427+0.000083. 
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Fig. 5. — Same plots as Figure H] except the limit of ellipticity for using objects, and the 
limit of this figure is \5\ <1. KSB method has a number density 2.48(/min 2 ) and obtained 
fitting parameter 0.31087+0.000737. E-HOLICs method has a number density 3.41(/min 2 ) 
and obtained fitting parameter -0.00447+0.000104. 




Fig. 6. — Same plots as Figure H] except the limit of for using objects, and the limit of 
this figure is 2.2< r^. KSB method has a number density 3.85(/min 2 ) and obtained fitting 
parameter 0.02777+0.000231. E-HOLICs method has a number density 4.05(/min 2 ) and 
obtained fitting parameter -0.00757-0.000233. 
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Fig. 7. — Same plots as Figure @] except the limit of SN for using objects, and the limit 
of this figure is 2<SN. KSB method has a number density 4.99(/min 2 ) and obtained fitting 
parameter -0.05397-0.000105. E-HOLICs method has a number density 5.38(/min 2 ) and 
obtained fitting parameter -0.06047-0.000553. 




Fig. 8.— Same plots as Figure 0] except the limits for using objects, and the limits of 
this figure is \5\ <1, 2.2< and 4<SN. KSB method has a number density 23.70(/min 2 ) 
and obtained fitting parameter 0.00997-0.000261. E-HOLICs method has a number density 
24.42(/min 2 ) and obtained fitting parameter -0.00147-0.000863. 
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